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We present a study on the modification of the electronic structure and hole-doping effect for 
the layered dichalcogenide WSe 2 with a multi-valley band structure, where Ta is doped on the 
W site along with a partial substitution of Te for its lighter counterpart Se. By means of band- 
structure calculations and specific-heat measurements, the introduction of Te is theoretically and 
experimentally found to change the electronic states in WSe 2 - While in WSe 2 the valence-band 
maximum is located at the F point, the introduction of Te raises the bands at the K point with 
respect to the T point. In addition, thermal-transport measurements reveal a smaller thermal 
conductivity at room temperature of Wi-^Ta^Sei.eTeo^ than reported for Wi_ :E Ta ;E Se 2 . However, 
when approaching 900 K, the thermal conductivities of both systems converge while the resistivity 
in Wi_a,Ta^Sei. 6 Teo .4 is larger than in Wi_ I! Ta a; Se 2 , leading to comparable but slightly smaller 
values of the figure of merit in Wi-iTaxSei.gTeo^- 

PACS numbers: 72.20.-i,72.20.Pa,71.20.-b 


I. INTRODUCTION 

Bulk dichalcogenides TX 2 with transition metals T 
and X = S, Se, or Te belong to the class of layered 
semimetals / semiconductors 1 which attracted consid¬ 
erable attention for decades. They exhibit a rich va¬ 
riety of physical propert ies, such as superconductivity 
in, e.g., NbSe 2 or TaS 2 p® coexisting but competing 
charge-density wave orderJMH an( f thermoelectrical de¬ 
vice functions.^ In many cases, these materials possess 
Fermi surfaces consisting of many valleys, which are 
somehow related to, or responsible for, these phenom¬ 
ena. 

The discovery of graphene has further increased the 
interest in layered materials with valley degrees of free¬ 
dom. It was found that the physical properties of bulk 
layered materials can be significantly altered when thin¬ 
ning them down to atomically flat layers.®® Natu¬ 
rally, different layered systems came into the focus of 
research. Among them are the monolayered dichalco¬ 
genides which are often regarded as two-dimensional 
semiconducting or semimetallic analogs to graphene: 1 ® 15 ! 
The electronic band structure plays a crucial role also 
in the physics of these thin materials with impact on 
spintronics,®optoelectronics devices,® r ® or the emerg¬ 
ing field of valleytronics, e.g., in M 0 S 2 , MoSe 2 , WS 2 , and 
WSe2.^ 20 ^31 The strong spin-orbit interaction in some 
transition-metal dichalcogenides and the lack of inversion 
symmetry in their monolayer variants lead to valley-spin 
coupling, adding a new feature to their zoo of exotic 
properties.® ' 22 24 2 - Another example of their intrigu¬ 
ing nature is the electric-field-induced superconductiv¬ 
ity with an optimum T c exceeding 10 K in an electric 
double-layer transistor structure made from M 0 S 2 .®! 31 ! 
The band structures of the aforementioned hexagonal 
multilayer (bulk) and monolayer WS 2 , WSe 2 , and M 0 S 2 


differ as well. While the bulk systems exhibit indirect 
band gaps, the monolayered counterparts have direct 
band gaps at the inequivalent K and K' poin ts at the 
corners of the hexagonal Brillouin zone. 15 16 32-341 There¬ 
fore, to control the band structure is one of the most 
important issues in the study of those TX 2 compounds. 

The presence of many valleys in the band structure is 
in general interesting also for thermoelectrical applica¬ 
tions since it is well known to enhance the thermoelectri¬ 
cal performance.® Thermoelectrical materials are also in 
the focus of current research because they offer the pos¬ 
sibility to transform thermal waste heat back into usable 
electrical power!®®!® A measure of the thermoelectric ef¬ 
ficiency of such materials is the figure of merit (FOM) 
ZT = S 2 T/(pk). It consists of the thermopower or See- 
beck coefficient S, the longitudinal resistivity p , and the 
total thermal conductivity n = n e \ + /c p h, where n e \ and 
Kph are the contributions of the mobile charge carriers 
and the lattice, respectively. The term S 2 /p is often re¬ 
ferred to as the power factor. To maximize ZT, a large 
thermopower and a small resistivity along with a small 
thermal conductivity are required. Conventional metals 
are usually not good candidates for high-efficiency ther¬ 
moelectrics since the Seebeck coefficient is generally small 
due to the Fermi degeneracy. Instead, doped semiconduc¬ 
tors are promising materials where it is possible to con¬ 
trol the charge-carrier concentration and hence the elec¬ 
trical conductivity and the Seebeck coefficient by doping. 
However, most of such systems suffer from a good lattice 
heat conductivity. The guiding principle canbe abbre¬ 
viated as “phonon glass + electron crystal”,® 38 i.e., a 
system which ideally consists of independent charge- and 
heat-transport channels to obtain low p and n values at 
the same time. In this context, bulk dichalcogenides at¬ 
tracted much interest when it was found that several of 
them exhibit at room temperature a small heat conduc- 
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tivity of only ~ 2 W/K m, among them \VSe2P Back in 
the 1960s and 1970s, this material was intensively studied 
by doping into both the W (cation) and the Se (anion) 
lattice sites®®! as well as by intercalation of metal ele¬ 
ments to bridge the chalcogen layers.®! 

WSe2 is a p-type semiconductor and crystallizes in the 
hexagonal PQ^/mmc structur e (sp ace group 194), usu¬ 
ally abbreviated as 2H- \VSe2 .EES It consists of trilayer 
building blocks Se-W-Se with strong covalent bonds. 
These blocks are separated by only weakly bonded van 
der Waals gaps. Each W ion is coordinated by six Se ions 
with a trigonal prism configuration. The unit cell con¬ 
sists of two formula units along the c axis. Several band- 
structure calculations can be found in literature, the most 
recent one in Ref. 1251 The valence-band maximum lies at 
the T point (zone center) almost degenerate with the only 
slightly lower-lying band at the K (and K') point. While 
the band dispersion at the T point has an almost isotropic 
nature, at the K point it is anisotropic. By doping pen- 
tavalent Ta 5+ on the W 4+ site, holes are introduced into 
the valence band, and the large electrical resistivity of 
pure WSe 2 is successfully suppressed: Wi_ a; Ta 2 ,Se2 be¬ 
comes metallic at small Ta concentrations x ~ 0.03 with 
room-temperature resistivities of the order ~ mflcm. Up 
to x ~ 0.35, the struct ure rem ains hexagonal PQ^/mmc 
with p-type conduction ! 9 ! 41 ! 43 ] On the other hand, there 
are fewer works published on WSe2- y Te, / , i.e., the sub¬ 
stitution of isovalent Te 2 ~ for Se 2- . The end mem¬ 
ber WTe2, crystallizing in an orthorhombic structure 
(Pmmm ), is a senrinretal with a negative Hall coeffi¬ 
cient. It seems that at least up to y = 0.5 the hexag¬ 
onal WSe2 structure is retainedP-EH As for the prepa¬ 
ration of these dichalcogenides, it should be noted that 
they do not melt congruently and dissociate at elevated 
temperatures. This was discussed earlier as a problem 
in achieving samples with high packing densities.® 4 Also, 
the physical properties were not very reproduciblef 41 ! 

Here we present a comprehensive study on 
Wi-xTaaSei.eTecu with 0 < x < 0.06 by means 
of band-structure calculations, transport, and thermo¬ 
dynamical measurements to elucidate the effect of Te 
doping on the electronic structure. Interestingly, we 
found evidence that the valence-band maximum shifts 
from the T point in WSe2 toward the I\ point when 
substituting Se by Te while keeping the crystal structure, 
which is reminiscent of the aforementioned situation 
in the monolayer dichalcogenides. The introduction 
of Te was also found to further suppress the thermal 
conductivity around room temperature. 


The paper is organized as follows: In the next sec¬ 
tion (Sec. 0 . the sample preparation and experimental 
and theoretical calculation methods are described . Afte r 
overviewing the electronic band structure (Sec. Ill A I, 
the basic transport properties and structural data be¬ 
low room temperature for Wi-ajTa^Sei.eTeo^ are intro¬ 
duced, to gethe r with those for W i_ x Taa;Se2 for compari¬ 
son (Sec. Ill B). Next (Sec. Ill C[ ) we focus on the change 
in the electronic states of WSe2 due to the Se replacement 


with Te by discussing the results of specific-heat data in 
the light of the band-structure calculations. In the latter 
half of this paper (Sec. HID), we present the thermo¬ 
electric properties of Wi-^Ta^Sei.eTeo^, i.e., thermal- 
transport data up to ~ 850 K from which the FOM is 
calculated. Section [TV] is devoted to the summary of the 
present work. 


II. SAMPLE PREPARATION AND METHODS 

Poly crystals!®! of Wi_! Ta^Sei. 6Teo. 4 for x = 0, 0.02, 
0.025, 0.03, 0.035, 0.0375, 0.04, and 0.06 were prepared 
in three steps. First, stoichiometric amounts of the ele¬ 
ments W (purity 99.9%), Ta (99.9%), Se (99.999%), and 
Te (99.999%) were thoroughly ground, mixed, sealed into 
evacuated quartz tubes, and kept for 48 h at 700°C. The 
resulting reaction product was free-flowing blackish pow¬ 
der. Second, the powder was reground, pelletized, again 
sealed into evacuated quartz tubes, and kept for 72 h to 
96 h at 1000°C. Third, the reaction product was ground 
again. Approximately 600 mg of the powder was used to 
synthesize a final batch by employing a high-pressure (p ) 
/ high-temperature (T) technique using a cubic anvil cell. 
The latter approach was chosen to overcome the afore¬ 
mentioned problem to achieve high-packing densities.®! 
Practically, the powder was first pressed at 2 GPa at 
room temperature. Next, the temperature was increased 
to about 1100°C. This temperature was kept for about 
10 minutes, then the material was thermally quenched. 
The pressure was released after the temperature had 
dropped back to room temperature. For each of the 
three steps, the reaction product was checked by x-ray 
diffraction (XRD). The targeted compounds had already 
formed after the first step although the respective XRD 
peaks were very broad. This probably indicates a disor¬ 
dered stacking of the characteristic (Se,Te) - W-(Se,Te) 
trilayers. The second and third reaction steps respec¬ 
tively lead to a sharpening of the XRD patterns. For the 
various measurements, rectangularly- and cylindrically- 
shaped samples were cut from these batches.® We note 
that we did not observe any cleavage-like behavior when 
cutting the samples. This is probably a result of the ap¬ 
plication of hydrostatic pressure in our cubic anvil press, 
leading to more isotropic samples without preferred ori¬ 
entation in spite of the layered structure of WSe2- As 
a control experiment, we also prepared pristine WSe2 
and Wi_ x Ta x Se 2 {x = 0.02, 0.03, 0.04, 0.05, 0.06) by 
the same high-p / high-T synthesis but without any pre¬ 
reaction, yielding sharp XRD line patterns. Compara¬ 
bly sharp line patterns for Wi_ x Taa,Sei.6Teo.4 were only 
achieved after the three-step growth process. We note 
that some batches contained small amounts of unreacted 
Se / Te after the high-p / high-T synthesis step. From 
TG-DTA (Rigaku ThermoPlus Evo TG 8120) analyses, 
we found that the samples lost some weight when heat¬ 
ing them up. At higher temperatures the system starts 
to dissociate in agreement with an old report.®! There- 
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fore, we restricted our high-temperature experiments to 
below approximately 850 K. Moreover, most of the ex¬ 
amined specimens were annealed 24 h at 150°C before 
the high-T measurements. 

Below 300 K, the longitudinal resistivity p xx and the 
Hall resistivity p yx were measured by a standard five- 
probe technique using a commercial system (Quantum 
Design, PPMS). Low-T longitudinal Seebeck S xx and 
thermal conductivity k xx datgP 5 were taken with two 
separate home-built setups, each mounted on a PPMS 
cryostat. Above room temperature up to approximately 
850 K, p xx and S were measured simultaneously in a 
ZEM-3 system (ULVAC Technologies), where the sam¬ 
ple is held in He atmosphere by two Pt or Ni electrode 
stamps acting as current leads. Two thermocouples were 
pressed to one sample surface acting as voltage pads and 
were used for the Seebeck measurement. High-T resis¬ 
tivity and thermopower data were taken upon cooling. 
The respective low-T measurements were carried out us¬ 
ing the same samples after the high-T experiment. The 
thermal conductivity above room temperature was esti¬ 
mated using the formula n = D c p d 300K with the ther¬ 
mal diffusivity D 1 the room-temperature density of the 
respective samples d 300K , and the specific heat c p . The 
thermal diffusivity was measured by employing the laser- 
flash method in a commercial Netzsch LFA-457 appara¬ 
tus. The density of the samples was estimated from their 
mass and dimensions. Specific-heat data below 300 K 
were measured by a relaxation-time method using the 
PPMS. The addenda heat capacity was measured before 
mounting the sample and eventually subtracted from the 
total signal. Since the specific heat of Wi_ x Ta x Sei. 6 Teo .4 
at room temperature already exceeds 95% of the classi¬ 
cal Dulong-Petit limit c DP , c p at higher temperatures was 
calculated as c p = yT + c DP with the Sommerfeld param¬ 
eter 7 , which was determined from c p vs T plots at low 

temperatures.^ 

The first-principles band-structure calculations were 
performed with the WIEN2k code employing the full- 
potential linearized augmented plane-wave method.^ We 
used the Perdew-Burke-Ernzerhof exchange-correlation 
functional.^! In the calculations for WSe 2 , the measured 
lattice parameters (a = 3.289 A, c = 12.988 A) and 
the positional parameter z = 0.129 (W-Se layer dis¬ 
tance) reported in Ref. [55] were employed. For hypo¬ 
thetical 2 H-WTe 2 , we optimized the parameters by us¬ 
ing the scalar-relativistic approximation'^! (a = 3.555 A, 
c = 14.447 A, and z = 0.126). Spin-orbit coupling 
was also taken into account as a relativistic effect in 
the band-structure calculations.^ The structure opti¬ 
mization (band-structure calculation) was carried out 
with the cutoff i?A' max = 8.5 (10.0) and 6x6x4 
(12 x 12 x 8 ) k points. The density of states was cal¬ 
culated with 36 x 36 x 8 k points using the tetrahedron 
method.^ 



FIG. 1: (Color online) Band structure calculated for (a), (c) 
2H-WSe2 and (b), (d) hypothetical 2H-WTe2- The values of 
the lattice constants used for the calculations are indicated 
in panels (a) and (b). For the details, see text. Panels (c) 
and (d) provide expanded views of the band dispersions along 
the k z directions from the F and K points of the hexagonal 
Brillouin zone. In each panel, the Fermi level is indicated by 
a horizontal line {E = 0 eV). 


III. RESULTS AND DISCUSSION 

A. Band Structure 

Figure [l] summarizes the results of band-structure cal¬ 
culations. In Fig. [I] (a), the band structure based on ex¬ 
perimentally determined lattice constants for 2 H-WSe 2 is 
plotted. The valence-band maximum is located at the T 
point which has an almost isotropic dispersion. 59 The sec¬ 
ond, almost degenerate valence-band maximum is found 
at the K and K' points and lies approximately 100 meV 
lower in energy, yielding a unique band structure as dis¬ 
cussed in Ref. [25] Panel (c) gives an expanded view of 
the band dispersion along the k z direction at the T and K 
points, i.e., T-A and K-H , respectively. The in-plane and 
out-of-plane band widths around the K and K’ points are 
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much more different than around the T point. Hence, the 
band dispersion around K or K' has an anisotropic, two- 
dimensional nature. 

To investigate how the WSe 2 band structure changes 
upon Te doping, we calculated the energy dispersion 
for hypothetical 2 H-WTe 2 based on optimized lattice 
and positional parameters. It is important to mention 
that hexagonal 2H-WTe2 does not exist in nature. In 
reality, WTe 2 crystallizes in the orthorhombic Pmmm 
structure! 60 * 61 * The result is shown in Fig. [l] (b). Ex¬ 
panded views along the k z directions from the T and K 
points are plotted in panel (d). Interestingly, the sub¬ 
stitution of the heavier Te for Se shifts the valence-band 
maximum from T to K (and K') 1 but the energy dif¬ 
ference between them remains small. Most importantly, 
this shift has implications on the density of states (DOS) 
at the onset of the hole band. While the band-edge DOS 
in 2H-WSe2 is dominated by the almost isotropic band at 
the T point, it originates from the more anisotropic bands 
at the K and K' points in hypothetical 2H-WTe2- As 
long as the 2H-type hexagonal structure is retained, it is 
reasonably expected that the relative position of the band 
maximum gradually and continuously changes upon in¬ 
creasing the Te concentration in WSe 2 _j / Te ?/ . The region 
dominating the band-edge DOS changes from around the 
T point to around the K and K' points, when the energy 
level at the local maxima for the T and K points be¬ 
come the same. Following our calculations this happens 
at y ~ 0.74. 


B. Sample Characterization 

Figures [2] (a) and (b) show the Ta-doping depen¬ 
dence of the lattice constants a and c, respectively, for 
Wi_ x Ta a ,Sei. 6 Te 0 . 4 ( y = 0.4; filled red symbols) and 
Wi_ x Ta x Se 2 (y = 0; open blue symbols). According to 
the introduction of Ta and Te in the two different lattice 
sites, the lattice constants change systematically with x 
and y 1 indicating the successful formation of solid solu¬ 
tions. Ta doping increases the a-axis length whereas the 
c axis shrinks. The increase in the absolute values be¬ 
tween y = 0 and y = 0.4 reflects the replacement of 20% 
of the Se 2- ions with larger Te 2- ions. Both of the Te 
and Ta ions readily substitute their counterparts. 

The charge-carrier concentration as estimated from 
Hall resistivity p yx measurements at 300 K increases from 
~ 8.7 x 10 19 cm- 3 for x = 0.02 to ~ 9.4 x 10 2 ° cm" 3 for 
x = 0.06. These data are plotted with filled symbols as 
charge-carrier concentration per W site n in Fig. [2] (c); 
n = 0.01 corresponds to a charge-carrier concentration 
of ~ 1.65 x 10 2 ° cm -3 . We also estimated n at 5 K 
(not shown), finding only small changes from the values 
at 300 K. At both temperatures, the Hall resistivity p yx 
is linear in the magnetic field, indicating that there are 
only hole-type charge carriers. We are not able to mea¬ 
sure the Hall resistivity at elevated temperatures, but 
for low-doped Wi_ x Ta a ;Se 2 , a temperature-independent 



x x 



FIG. 2: (Color online) (a), (b) Evolution of the lattice con¬ 
stants a and c with x for Wi_ a! Ta ;E Se 2 (open blue symbols; 
y — 0) and Wi- t Ta^Sei.6Teo .4 (closed red symbols; y = 0.4). 
(c) Actual charge carrier concentration per W site n as es¬ 
timated from Hall resistivity p yx at 300 K is plotted against 
the nominal Ta-doping concentration x of Wi_a,Ta x Sei. 6 Teo .4 
(filled symbols). Data for Wi_ :c Ta x Se 2 (open symbols) are 
included for comparison; see text. The dotted line indicates 
n = x. The color code of the different data points is the same 
as in Fig. [4] 


charge-carrier concentration up to T ~ 900 K was re¬ 
ported earlier. 43 For comparison, respective data for 
Wi_ x Ta x Se 2 are also shown in open symbols. While 
the actual charge-carrier concentrations n are close to 
the nominal doping levels x for Wi_ x Ta a ;Se 2 , this is not 
the case for x < 0.04 in the Wi_ x Ta x Sei. 6 Te 0 .4 system. 
We find a systematic suppression of n compared to the 
nominal carrier concentration corresponding to the Ta- 
doping level x, i.e., one hole per Ta 5+ . Possibly, some of 
the introduced charge carriers are annihilated in the low- 
doped samples due to compensation effects arising from 
deficiency in the Se - Te stoichiometry. Such an effect is 
absent in the “pure” diselenide system. This annihilation 
does not seem to be active when a certain amount of Ta 
is doped and the system becomes metallic. In Ref. RE5I 
the author reports a smaller than nominal charge-carrier 
concentration also for Wi_ x Ta x Se 2 and speculates about 
a compensation effect due to impurities at low doping 
concentrations. 

Figure [3] summarizes the resistivity p xx data below 
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FIG. 3: (Color online) Temperature dependence of the resis¬ 
tivity p xx below room temperature of Wi_ x Ta ;E Se 2 (dashed 
blue lines; y = 0) and Wi- t Ta^Sei. 6 Teo .4 (solid red lines; 
y = 0.4) for 0 < x < 0.06. Note the change in the scale of the 
vertical axis. 


300 K for both the S e-unsubstituted Wi_ a; Ta a ;Se 2 (dashed 
blue lines) and the Se -substituted Wi_jTa^Sei. 6 Teo .4 
(solid red lines). In both systems, the partial replace¬ 
ment of W 4+ by Ta 5+ introduces holes and systemati¬ 
cally lowers the resistivity with x. Already a small Ta 
concentration of x < 0.02 drastically su ppres ses p xx by 
orders of magnitude compared to x = 0 . 4 d®l WSe 2 and 
WSei. 6 Teo .4 have room-temperature resistivities of the 
order of ~ 10 flcm whereas p at 300 K of Wo.ggTao.oiSe 2 
has dropped to ~ 0.1 flcm, and Wo.gsTao.o 2 Se 2 exhibits 
already less than 10 -2 flcm. For Wi_ K Ta x Sei. 6 Teo. 4 , 
the resistivity drops for x = 0.03 below 10 -2 flcm. 
Metallic behavior is observed in both series above x > 
0.04, although the absolute values of p xx remain in 
the mflcrn range. The resistivity of the doping series 
Wi_ a ,Ta a ,Se 2 (see also Refs. HD and H51) drops faster than 
in Wi_ x Taa;Sei. 6 Teo. 4 . This is probably due to addi¬ 
tional disorder caused by the partial substitution of Te 
for Se. We note that in Ref. [4D especially for x = 0.01, 
smaller values are reported for p xx , whereas those of the 
metallic compositions are comparable. 



FIG. 4: (Color online) Specific-heat c v data of 

Wi_a: Ta^Sei. 6 Teo .4 (filled symbols) and Wi- a: Ta a; Se 2 
(open symbols): (a) c v vs T as measured. The dashed line 
depicts the specific heat calculated in the Debye model plus 
electronic contribution for x = 0.06; see text. The inset 
shows an expanded view around 100 K displayed as c p /T vs 
T. (b) The same data plotted as c p /T vs T 2 at low T. The 
lines therein are linear fits to the data. An expanded view of 
the extrapolation to 0 K is shown in (c). Solid lines refer to 
Wi_a:Ta^Sei. 6 Teo. 4 (y = 0.4), dashed lines to Wi_ I ,Ta a; Se 2 
(y — 0). The intercepts at T = 0 K give the electronic 
specific-heat coefficients 7 ; see text. 


C. Specific Heat 


Next, we discuss the Te-doping effect in terms of 
specific-heat data below T ~ 330 K down to approx¬ 
imately 1.9 K. There is no specific-heat study on the 
doped systems reported in the literature. For the mother 
compound WSe 2 , we found only one publication report¬ 
ing c p data above 60 Here we show that the analy¬ 
sis of the Cp data provides important information about 
the change in the electronic states when going from 
Wi_ a ,Ta a; Se 2 to Wi-a,Ta^Sei.@Teo. 4 . These data are 
summarized in Fig. [4] for Wi-^,Ta^Sei. 6 Te 0 .4 (x = 0.02, 
0.03, 0.04, and 0.06; filled symbols) along with data for 
Wi_ a; Ta a ;Se 2 (x = 0, 0.02, and 0.06; open symbols). In 
Fig. [4] (a), the data are displayed as c p vs. T. The specific 
heat of all samples is very similar on this scale. Around 
room temperature, each sample has released more than 
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~ 95% of the entropy expected in the classical Dulong- 
Petit limit c DP plus its respective electronic contribution 
7 T. This is exenrplarily indicated as a dashed line in 
Fig. [4] (a) for x = 0.06. 

However, the expanded plot shown as c p /T vs T in 
the inset of Fig. [I] (a) reveals that there is a clear impact 
on c p when changing x and y. The difference between 
c p of Wi_ ;c Ta ;c Se 2 (y = 0 ) and Wi_ x Ta x Sei. 6 Te 0 .4 (y = 
0.4) is due to the change in the phonons (lattice) caused 
by the replacement of 20% Se with Te. Moreover, we 
could also successfully resolve very small but systematic 
changes in both systems due to the increase of the Ta 
concentration x, and hence a change of the electronic 
specific heat. This can be seen in Fig. |4] (b), which shows 
the specific heat at low temperatures displayed as c p /T 
vs T 2 , where the phononic contribution is small and the 
electronic sector can be studied. Figure[4](c) provides an 
expanded view of the extrapolation of the low-T data to 
0 K (linear fits to the respective data). Solid lines refer 
to Wi_a;Taa;Sei. 6 Teo. 4 , dashed lines to Wi_ a; Ta 2 ,Se 2 . 

Conventional Debye fits to the data c p vs T for T <5 K 
using 

c p = c e i + c ph = 7 T + A 3 T 3 (1) 

yield good descriptions with the electronic specific-heat 
coefficient 7 and the coefficient of the phononic part A 3 . 
From the latter, the Debye temperatures ©d of each sam¬ 
ple were calculated via A 3 = (12/5) 7 r 4 ./V N a&b /©d with 
the number of atoms per formula unit N = 3, the Avo- 
gadro number JVa, and the Boltzmann constant fee- For 
all measured samples of the Wi-^Ta^Sei.gTeo^ system, 
we found ©d ~ 224 K and for the Wi_ a ,Ta a; Se 2 system, 
©d ~ 235 K for x = 0, 0.02, 0.06, and 211 K (x = 0.01), 
and 228 K (x = 0.04). The difference in the released 
entropy when going from the pristine diselenide system 
( y — 0) to the Se-substituted system (y = 0.4) up to 
approximately 200 K where the different c p data start to 
converge amounts to approximately 3 J/(mol K). 

The undoped compounds without Ta exhibit an one- 
order-of-magnitude smaller electronic contribution to 
the specific heat as already suggested by the insulat¬ 
ing nature of these specimens with a band gap. For 
both systems we observe a systematic change of 7 = 
( 7 r 2 fcg/ 3 )DOS with the actual charge-carrier concentra¬ 
tion n as plotted in Fig. [5] (a). This reflects the change 
in the DOS upon hole doping, clearly showing that the 
filling dependence of 7 appears to be qualitatively dif¬ 
ferent between the two systems. In view of our band- 
structure calculations shown in Fig. [I] the substitution 
of Te for Se leads to a rise of the bands at the K point 
with respect to the T point. Since the band dispersion 
at the K point is highly anisotropic whereas it is almost 
isotropic at the T point, charge carriers in Wi_ x Ta a; Se 2 
and Wi_ x Taa;Sei. 6 Teo .4 upon Ta doping may be accom¬ 
modated into bands of different degree of anisotropy. Mo¬ 
tivated by this, we tried to fit the 7 ( 71 ) data points to 
different formulas: (i) for Wi_ a ,Ta a; Se 2 to 

7 ( 74 )= an 1 / 3 ( 2 ) 



FIG. 5: (Color online) (a) Charge-carrier concentration de¬ 
pendence of 7 of Wi—xTax Sei .6 Teo .4 (filled symbols) and 
Wi_ x Taa,Se 2 (open symbols). The dashed (red) and dashed- 
dotted (blue) lines are fits to the data assuming differences 
between the two systems in the underlying band structure; see 
text. The color code of the different data points is the same as 
in Fig. [I] Since the data points for both mother compounds 
without Ta are overlapping on this scale, the one referring to 
WSe 2 is plotted with a smaller symbol, (b) Calculated density 
of states of Wi_ I! Ta a; Se 2 and hypothetical 2 H-Wi_a;Taa:Te 2 
as a function of charge-carrier concentration per W site n. 
The arrow marks the onset of filling of the isotropic band at 
the T point for hypothetical 2 H-Wi_ a: Taa;Te 2 ; see text. 

as expected for isotropic bands, and (ii) for 
Wi-xTa^Sei.eTeo^ to 

7 (n) = 70 + bn 1/3 . (3) 

Here, 70 is an offset which accounts for the constant 
DOS of an ideally two-dimensional band structure. 
These approaches yield the dashed-dotted (Wi_ a; Ta a ;Se 2 ) 
and dashed (Wi_ 2 ,Ta^Sei. 6 Teo. 4 ) lines in Fig. i (a), 
well accounting for a qualitative difference between 
both systems. For Wi-a, Ta^Sei.eTeo. 4 , we find 70 = 
0.52 mJ/molK 2 . We note that for Wi_ 2 ; Ta a ;Se 2 the 7 
value of the undoped compound WSe 2 was included to 
the fit without problems whereas for Wi-^Ta^Sei.gTeo^, 
data for x = 0 had to be excluded to obtain good fit¬ 
ting results, which is another indication for the differ¬ 
ence in the DOS of both systems. Since the bands 
at the r and K points change their relative position 
when going from WSe 2 to WTe 2 , our observation is 
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reasonable. Upon increasing the Te concentration, the 
bands at both points are getting closer in energy: In 
Wi_ x Ta x Sei. 6 Te 0 . 4 , charge carriers also populate the 
bands of highly anisotropic character at the K point. 
Hence the constant offset y 0 is needed to describe 7 (n) 
properly. 

To further strengthen this argument, the calculated 
DOS as a function of charge-carrier concentration per W 
site n is shown in Fig. [5] (b). Again the dashed-dotted line 
refers to Wi_ x Ta x Se 2 and the dashed line to hypothetical 
2H-Wi_ x Ta x Te 2 . The ditelluride system exhibits an ini¬ 
tial steplike increase of the DOS as expected for the filling 
of charge carriers into a highly anisotropic band at the 
K point. The arrow in Fig. [5] (b) marks the onset of fill¬ 
ing into the isotropic bands at the T point. By contrast, 
the diselcnide system exhibits a smoother initial filling, 
indicating the isotropic character of the band structure 
at the T point. A direct comparison of Figs. [5] (a) and 
(b) is difficult, since Wi_ x Ta x Sei. 6 Teo .4 contains Se and 
Te and hence is located in between the two extreme cases 
shown in panel (b). However, the filling dependence of 7 
for Wi_ x Ta x Sei. 6 Te 0 .4 can be well accounted for if the 
onset of the filling into the isotropic band at the T point 
(as indicated by the arrow) nearly coincides with n = 0. 


D. Thermoelectric Properties 

Having confirmed that the band structure is suc¬ 
cessfully modified by the partial replacement of Se, 
next we studied the thermoelectric properties of the 
Wi_ x Ta x Sei. 6 Teo .4 system. In Fig. [6j the tempera¬ 
ture dependence of transport data up to ~ 850 K for 
Wi_ x Ta x Sei. 6 Teo .4 is summarized: (a) resistivity p xx , 
(b) thermopower S, and (c) the corresponding power fac¬ 
tor S 2 / p xx . Although different experimental setups were 
used for the low-T and high-T measurements, respective 
resistivity data sets agree with each other around room 
temperature within the experimental error bars, except 
for x = 0.04. For x < 0.04, each sample exhibits a 
minimum in p xx at elevated temperatures. Above about 
650 K, all of them feature a metal-like positive slope of 
Pxx against T. The low-T and high-T Seebeck coeffi¬ 
cients shown in panel (b) also agree well around 300 K. 
The Ta doping leads to a very systematic change of S , 
too. For all x , the thermopower increases with T up to 
the highest measurement temperature 850 K, although 
the slope flattens above room temperature for most of 
them. For degenerated semiconductors, the thermopower 
is expected to be linear in temperature, which we in¬ 
deed observe below the flattening. The most metallic 
sample x = 0.06 retains the linearity up to the highest 
measurement temperature. At around 850 K, the largest 
thermopower S ~ 300 /xV/K is found for the low-doped 
samples. From the data shown in Fig. [6] (a) and (b), the 
power factors S 2 /p xx were calculated for each sample and 
plotted in (c). We observe a maximum in the power fac¬ 
tor of S 2 /p xx ~ 8 /xW/K 2 cm at around 850 K for samples 



FIG. 6: (Color online) Transport data of Wi- x Ta x Sei. 6 Teo ,4 
for 0.02 < x < 0.06. (a) Resistivity p xcc , (b) thermopower S, 
and (c) power factor S 2 / p xx are plotted against temperature 
up to 850 K. In (a), data for WSei.6Teo.4 is also included for 
comparison. 


in the intermediate doping range x = 0.03 — 0.04. This 
is due to the benefit of an increased conductivity while 
the thermopower is still fairly large even for the more 
metallic specimen with x = 0.04. 

The thermal conductivity and respective FOM ZT = 
S 2 T/(pn) for these samples are plotted in Figs. 0 (a) 
and (b), respectively. For the thermal-conductivity mea¬ 
surements, not only the experimental setups differ be¬ 
tween high-T and low-T measurements, but also differ¬ 
ent samples were used . 50 The agreement between low-T 
and high-T data is less satisfactory than in the resistiv¬ 
ity and thermopower measurements, but still within the 
acceptable range. Except for x = 0.03, all data exhibit a 
maximum in k(T ) below 100 K. Toward higher temper¬ 
atures, the thermal conductivity for all x continuously 
decreases due to the three-phonon-scattering process and 
















FIG. 7: (Color online) Temperature dependence of (a) ther¬ 
mal conductivity k and (b) dimensionless thermoelectric fig¬ 
ure of merit ZT of Wi_a;Ta^Sei.gTeo .4 below 850 K. 


does not show any signature of saturation even at 850 K. 
The corresponding FOM plotted in Fig. [7] (b) increases 
with increasing temperature for all x. Above 800 K a 
maximum ZT « 0.35 is found for x = 0.035. 

For a better comparison, we replotted various prop¬ 
erties of Wi_ x Ta x Sei. 6 Te 0 .4 against the actual charge- 
carrier concentration n at 300 K (filled symbols) and at 
800 K (open symbols) in Fig. (a) p xx , (b) Hall mo¬ 
bility p = (nepxx)^ 1 with the element charge e, (c) S, 
(d) S 2 /pxx, (e) k, and (f) ZT. The dotted and dashed 
lines are guides to the eyes. The absolute values of the 
resistivities are not so different at 300 K and 800 K. The 
mobilities are plotted in Fig. [8j)b). Since Hall-effect mea¬ 
surements were carried out only at 300 K and not at ele¬ 
vated temperatures, panel (b) contains only 300 K data. 
The suppression of p xx with x is almost fully ascribed 
to the increase in n: An almost constant mobility is ob¬ 
served in this system. The thermopower at 300 K and 
800 K given in Fig. [8] (c) decreases upon Ta doping, in 
agreement with the decrease in the resistivity. The dif¬ 
ference AS = S $00 k ~ <5300 k in the absolute values for 
each sample amounts to around 100/iV/K. The maxi¬ 
mum of the power factor is clearly seen in Fig. H(d) 
between n = 0.02 and 0.04 at both temperatures. The 
absolute values increase by a factor of three or four when 
going from 300 K to 800 K. At the same time the av¬ 
erage thermal conductivity drops roughly to half of its 
room-temperature value, as depicted by the dotted and 
dashed lines in Fig. [8] (e). This reduction of the ther¬ 
mal conductivity should be attributed to an increased 



n n 


FIG. 8: (Color online) Summary of transport data for the 
Wi_a;Ta^Sei.6Teo .4 system at T = 300 K (filled symbols) and 
800 K (open symbols), (a) Resistivity, (b) Hall mobility, (c) 
thermopower, (d) power factor, (e) thermal conductivity, and 
(f) FOM are plotted as a function of charge-carrier concen¬ 
tration per W site n. The lines in the panels are guides to the 
eyes. The color code of the different data points is the same 
as in Fig. [4] 

phonon-phonon scattering rate. The dependence of the 
FOM on n shown in Fig. [ 8 ] (f) resembles that of the power 
factor and also exhibits a maximum in the intermediate 
Ta-doping range. The increase in ZT between 300 K and 
800 K amounts to more than one order of magnitude. 

Finally, we compare the results on Wi_ x Ta^Sei. 6 Teo .4 
with those on Wi_ a; Ta :E Se 2 reported in the literature. 
There is one paper focusing comprehensively on the ther¬ 
moelectric properties in Wi_ a; Ta a ;Se 2 by BrixnerP^ Ad¬ 
ditional transport data from the same group are sum¬ 
marized in Ref. [43] In Ref. HU, the temperature de¬ 
pendences of p , S, k, and ZT are shown for x = 0.01 
and 0.03. The respective values for the FOM at 800 K 
are ZT ss 0.6 for x = 0.01 and ZT k, 0.35 for 
x = 0.03, exceeding the maximum ZT value of 0.32 re¬ 
ported here for Wi_ 2 ,Ta^Sei. 6 Te 0 . 4 . The thermoelectric 
parameters for these samples are summarized for com¬ 
parison in Table [T] There are several differences between 
these systems: (i) For small x , the resistivity is smaller 
and the thermopower is larger in Wi_ x Ta a; Se 2 than in 
Wi-xTffeSei.eTetM, leading to larger power factors, (ii) 
Due to the disorder-scattering-induced reduction of the 
phononic mean-free path, the thermal conductivity in the 
Se-substituted system with y = 0.4 is clearly smaller at 
room temperature than the value of Wi_ a; Ta: C Se 2 (y = 0 ) 
reported in Ref. [41] Apparently, this merit of the Te 
codoping is overcompensated by the increase in the resis¬ 
tivity due to additional disorder in the anion sublattice, 
which increases the scattering rate of the p-type charge 
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TABLE I: Thermoelectric properties of Wi_ rc Ta a; Se 2 (from Ref. |41|) and Wi-^TaxSei^Teo^ (this work). The actual charge- 
carrier concentration n is presented in holes per W site, the resistivity in mficm, the thermopower in ^iV/K, the power factor 
in /iW/K 2 cm, and the thermal conductivity in W/K in; see text for details. 






Wi_ 2 ,Ta a ,Se 2 




Wi. 

—x TcLe Sg 1 _ 6 Teo . 4 


n 

T 

P 

S 

S 2 /p 

fctot 

ZT 

P 

S 

S / p H tot 

ZT 

0.006 

300 K 

7.0 

112 

1.8 

6.87 

0.01 

23.0 

194 

1.6 4.11 

0.01 

800 K 

6.9 

335 

16.3 

2.56 

0.60 

15.2 

302 

6.0 2.00 

0.24 

0.027 

300 K 

2.0 

90 

4.0 

7.31 

0.05 

5.1 

105 

2.2 3.10 

0.02 

800 K 

3.1 

166 

8.9 

2.23 

0.35 

5.8 

210 

7.6 1.90 

0.32 


carriers, (iii) The thermal conductivity® at 800 K in 
Wi_ a; Ta 2 ,Se 2 is larger than 2 W/K m, but almost compa¬ 
rable to Wi_a,Ta^Sei.gTeo. 4 , and hence at high tempera¬ 
tures the FOM is larger in the “pure” diselenide system. 

To see whether the phononic contributions to the total 
thermal conductivity are different, we assume that the 
Wiedemann-Franz law holds in these systems, and cal¬ 
culate the contribution of the charge carriers via K e i = 
Lq(T/ p) with the Lorenz number Lq = 2.44x 10 -8 V 2 /K 2 
of the Drude-Sommerfeld model. By subtracting n e \ from 
the total thermal conductivity Ktoto the phononic part 
Kph was calculated. In Ta^Sei.gTeo. 4 , the elec¬ 

tronic thermal conductivity n e i increases linearly with 
n (not shown) and amounts at room temperature to 
~ 0.03-0.3 W/K m and at 800 K to ~ 0.1-0.5 W/K m. 
This is less than 22% of the total thermal conductivity, 
and most of the thermal conductivity is attributed to the 
phononic contribution. From the data given in Ref. I4R 
we estimated the respective electronic and phononic ther¬ 
mal conductivities in Wi_ a; Ta a; Se 2 . To readily compare 
these numbers to Wi-xTa^Sei^Teo^, one has to take 
into account the real charge-carrier concentration per 
W site as estimated from Hall-effect measurements. In 
Ref. |TT]the hole concentration is only given for x = 0.01 
which corresponds to 0.004 per W site. Following Ref.PSJl 
the carrier concentration for x = 0.03 is 0.02 per W 
site. Using these numbers, the electronic and phononic 
contributions of the thermal conductivities from Ref. HU 


TABLE II: Electronic K e i and phononic K p h contributions 
to the total thermal conductivities Ktot = K P h + K e i in 
Wi- x Ta a; Se 2 (from Ref. l4ll) and W 1-1 Ta^Sei. 6 Teo .4 (this 
work). The actual charge-carrier concentration n is presented 
in holes per W site and the thermal conductivity values in 
W/K m; see text for details. 




Wi- 

— x T&a; SG 2 

Wi- 

— x TcLx Sg 1 .6 Tgo . 4 

n 

T 

^-el 

^ph 

K tot 

Kel 

ftph 

K'tot 

0.006 

300 K 

0.10 

6.77 

6.87 

0.03 

4.08 

4.11 

800 K 

0.28 

2.28 

2.56 

0.13 

1.87 

2.00 

0.027 

300 K 

0.36 

6.95 

7.31 

0.14 

2.96 

3.10 

800 K 

0.64 

1.59 

2.23 

0.33 

1.57 

1.90 


have to be compared with the respective values of our 
Wi-^Ta^Sei.gTeo.i samples for nominal x = 0.02 (n « 
0.006) and x = 0.035 (n « 0.027); see Fig. [ 2 ] (c). Ta¬ 
ble m summarizes the different contributions for these 
two charge-carrier concentrations at 300 K and 800 K. In 
both systems, the thermal conductivity is dominated by 
the phononic contribution. However, the phononic ther¬ 
mal conductivity is smaller in Wi_jTa^Sei. 6 Teo .4 than 
in Wi_ a; Ta a ;Se 2 , especially at 300 K. Hence both com¬ 
ponents of the thermal conductivity n = K p h + K e i are 
substantially reduced in Wi_jTa^Sei.gTeo. 4 . It might be 
promising to further change the Te concentration and see 
whether one can optimize the thermoelectric properties 
and exceed the FOM reported for Wi_ x Ta a; Se 2 . 


IV. SUMMARY 

We present a comprehensive theoretical and experi¬ 
mental study on Wi-j,Ta^Sei. 6 Teo .4 for small Ta concen¬ 
trations 0 < x < 0.06. From band-structure calculations 
and an analysis of the specific heat, we find clear evidence 
that upon introducing Te into the Se sites in WSe 2 , the 
band structure changes. The isotropic band at the T 
point is lowered in energy while the anisotropic bands at 
the K and K' points shift toward the Fermi level, leading 
to a change in the doped-hole character. This crossover 
was monitored in the electronic specific-heat coefficient 
which reflects the apparent change in the filling depen¬ 
dence of the density of states, namely, from filling charge 
carriers solely into the isotropic band at the T point in 
Wi_ a ,Ta a; Se 2 toward filling into both the isotropic band 
at the F point and the anisotropic bands at the K and 
K' points in Wi_!Ta^Sei. 6 Teo. 4 . 

Thermal- and electronic-transport measurements up 
to 850 K on this system yield that the maximum ther¬ 
moelectric figure of merit is about 0.3 at 850 K in the 
doping range 0.03 < x < 0.035. In comparison to 
the Wi_ a; Ta a: Se 2 system, Te doping was found to suc¬ 
cessfully suppress the thermal conductivity, especially 
around room temperature. However, this merit is over¬ 
compensated by an increased resistivity due to the ad¬ 
ditional disorder in the anion sublattice, which leads to 
stronger scattering of the hole carriers. 
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